Spin Waves in Random Spin Chains 
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We study quantum spin-1/2 Heisenberg ferromagnetic chains with dilute, random antiferromag- 
netic impurity bonds with modified spin-wave theory. By describing thermal excitations in the 
language of spin waves, we successfully observe a low-temperature Curie susceptibility due to for- 
mation of large spin clusters first predicted by the real-space renormalization-group approach, as well 
as a crossover to a pure ferromagnetic spin chain behavior at intermediate and high temperatures. 
We compare our results of the modified spin-wave theory to quantum Monte Carlo simulations. 



One-dimensional (ID) quantum Heisenberg sys- 
tems with random ferromagnetic-antiferromagnetic (FM- 
AFM) couplings are of interest due to the interplay be- 
tween quantum fluctuations ancLdisorder. An archetyp- 
ical system is SraCuPti-jjIrjjOeJj alloy of the pure com- 
pounds SrsCuPtOg (antiferromagnet) and SraCuIrOg 
(ferromagnet) . Another example is organic radical al- 
loy (p-CDpOV)i_a;(p-BDpOV)^.0 These systems are of- 
ten refered to as quasi-lD random spin chains since 
the interchain coupling cannot be neglected at low tem- 
peratures where three-dimensional (3D) ordering domi- 
nates.El Curie-like magnetic susceptibility with impurity- 
concentration dependent-.|Curie constant has been re- 
ported for these systems.EJj3 

Theoretically, random spin chains have beeii stud- 
ied by high-temperature series expansions Jflu quan- 
tum transfer matrix^|mf;thod,uiJ quantum Monte Carlo 



(QMC) simulatictfis,ErLl'l3 density matrix renormalization 
group methodp'E3 and most noticeably^ fl|e real-space 



renormalization-group (RSRG) methodO-B In a simple 
RSRG picture, the Curie-like susceptibility occurs be- 
cause ferromagnetic couplings lead to spin correlations 
that form spin clusters whose average size (N) and ef- 
fective spin (Seff) grow stochastically {N — AS^^j) at 
low temperatures. The coefficient A depends on the dis- 
tribution of random couplings, and, in general, increases 
with concentration of ferromagnetic couplings. There- 
fore, the magnetic susceptibility per spin (x) varies with 
temperature (T) as 
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i.e. Curie-like, in the low-T limit. Here, we assume /i — 
ks — 1 for convenience. 

In the RSRG treatment, one decimates the spins that 
are coupled by the strongest coupling in the system, and 
renormalizes the couplings among the remaining spins 
perturbatively. This procedure generates effective spins 
with effective (in general weaker) couplings, so that ex- 
citations at lower and lower energy scale can be probed. 
It becomes obvious that an initial broad distribution of 



spin couplings ensures the quick approach to the asymp- 
totic renormalization-group (RG) behavior. However, it 
is not the case for a Heisenberg spin chain with nearest- 
neighboring exchanges -I- J or — J with certain probabili- 
ties (a minimal model for reaLaystems) . In fact, an RSRG 
study of random spin chainsEHl has suggested the possibil- 
ity of having a crossover region spanning more than five 
orders of magnitude in T before the true RSRG scaling 
regime can be reached. This suggests that the crossover 
physics may be more relevant to experimental findings. 

In this paper, we model a random spin chain {S — 1/2) 
by the following Heisenberg Hamiltonian 
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where the coupling strength J,; is randomly chosen be- 
tween Jaf (AFM) with probibilityp and —Jp (FM) with 
probability 1 — p (see Fig. |l|). We are primarily inter- 
ested in the dilute doping regime, i.e. p <C 1, where RG- 
predicted Curie behavior is expected to occur at lower 
T with decreasing -p. In the absence of AFM couplings 
{.Jaf — 0), the system is a collection of independent FM 
spin segments with exponentially distributed numbers of 
spins. For small p, the static magnetic susceptibility x of 
the isolated spin segments shows a crossover from an or- 
dinary Curie behavior at high temperatures {T > Jp) to 
a characteristic x oc This FM-chain behaidor has 

been revealed by the modified spin- wave theory.tS Finite 
amount of Jaf introduces correlations between FM seg- 
ments, which form spin clusters that grow, in a random 
walk fashion, both in size and total spin at decreasing T. 
Statistical analysis shows that the RG-predicted Curie 
susceptibility 
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is expected at low T. Therefore, for dilute and not-too- 
weak AFM couplings, two crossover between distinctive 
T-dependence can be observed in this system. 
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FIG. 1: (a) FM segments with AFM couplings (double 
bonds). The segment length A''i, N2, are exponentially 
distributed with an average value (A'^) = 1/p. (b) Distribu- 
tion of FM and AFM couplings, (c) In the MSW theory, 
local chemical potentials (/^i, fi2, ■■■) are introduced for all 
FM segments, determined by zero-magnetization constraints 
in all segments self-consistently. 
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Unfortunately, by generating effective couplings per- 
turbatively, the RSRG method keeps discarding high- 
energy excitations and, thus, is not expected to describe 
correct physics well before the onset of the true RG be- 
havior. To address these higher energy scales (at in- 
termediate and high temperatures), we have general- 
ized the modified spin-wave theory (MSW) to random 
spin systems, with dominating ferromagnetic couplings, 
to study the crossover behavior from high-temperature, 
short-distance physics to low-temperature, long-distance 
physics. Our data suggests an unambiguous crossover 
from FM spin-chain to random spin-chain behavior, 
which becomes wider for smaller p. Comparing with 
our QMC sinmlations using continuous time loop clus- 
ter algorithm,E3 we discuss the limitations of the MSW 
method. 

The technical details of the application of the MSW 
theory-ill random spin chains have been presented else- 
where.E3 We summarize the MSW approach below. The 
main idea of the MSW theory is to introduce a chemi- 
cal potential (equivalent to a uniform magnetic field) for 
the spin-wave excitations to ensure zero-magnetization 
of a FM spin chain at finite temperatures, so that the 
number of spin waves no longer diverges in low dimen- 
sions.E3 The introduction of the chemical potential leads 
to a surprisingly good quantitative description of the 
thermodynamics of the FM spin chain at both low and 
high T. In the case of random spin chains, we intro- 
duce, as illustrated in Fig. |](c), one such chemical po- 
tential per FM segment, determined self-consistently by 
the zero-magnetization constraint of each individual FM 
segment. We apply a periodic boundary condition to 
each FM segment, so that zero-magnetization of a seg- 
ment implies zero-magnetization for every spin. With the 
physical significance of restoring the rotational symme- 
try broken during the introduction of spin waves, these 
zero-magnetization constraints ensure not only the finite 
number of spin waves at finite T, but the appropriate 
correlations between spin waves in different FM segments 
as well. A bosonic version of the generalized Bogoliubov 
transformation has been developed for this approach. 

Figure ^ (left panel) shows the static magnetic suscep- 
tibility per spin x for ^ random Heisenberg spin chain 



FIG. 2: Left panel shows static magnetic susceptibility per 
spin X of a 600-spin chain with 59 randomly distributed AFM 
couplings {Jaf = 0.5Jf). On the log-log scale, the two 
dashed lines with slope unity are exact results of the ordi- 
nary Curie susceptibility of independent spins at high T and 
the low-T RG-predicted Curie susceptibility (x = 0.75/r). 
The dotted line is the low-T exppflsion of x of a pure spin- 
1/2 FM chain (of infinite length) jHl approaching in the 

low-T limit (Eq. Right panel shows 7-sample averaged x 
for a shorter chain of 300 spin in 30 FM segments. High-T 
ordinary Curie behavior is not shown. 



of 600 spins, distributed in 30 FM segments of aver- 
age length 10 spins, i.e. p — 0.1, coupled by 59 AFM 
couplings with Jaf — 0.5 Jp. Three distinctive tem- 
perature regimes can be identified in the curve. For 
T / Jp > 3, X can be fit to an ordinary Curie law of in- 
dependent spins. (We neglect this regime in later plots.) 
For 0.2 < T I Jp < 3, x rises with decreasing T, foUowiag 
the trend of a infinitely-long, pure FM spin-1/2 spinjl3 
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This implies that spins start to correlate, forming inde- 
pendent FM segments. The AFM couplings are still too 
weak to affect thermodynamics of the random spin chain 
at T > 0.2Ji^. Below T = 0.2^^, x is in good agreement 
with the low-T, RG-predicted Curie behavior (Eq. || for 
p = 0.1). It is worth pointing out that we can identify the 
crossover temperature from the FM spin-chain physics to 
random spin-chain physics as the cross point of Eq. |^ and 
In the dilute doping limit (small p), we have 
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For p = 0.1, Tx = O.ITJf. The right panel in Fig. | 
shows averaged x for 7 random realizations of 300 spins 
in 30 segments with Jap — 0.5Jp. In this study, we 
require that the total magnetization of the corresponding 
classical ground state be zero, so x vanishes at low-T 
limit. This is responsible for the drop below T < 0.02, 
where there are no more segments to carry out RSRG 
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FIG. 3; Avaraged susceptibility x of random Heisenberg spin 
chains of 360 spins in 30 FM segments, with AFM impurity 
coupUngs Jaf = 0.5 Jf, for the MSW method (14 samples) 
and the QMC simulations (28 samples). The dashed line (x = 
ll/12r) is the low-T Curie susceptibility expected by the 
RSRG approach (Eq. ^). The inset shows the same data on a 
log-log scale. The dotted line is the low-T expansion of _x of 
a pure spin-1/2 FM chain (of infinite length)J13 i.e. Eq. 



to lower T. Nevertheless, the crossover in the vicinity of 
Tx — 0.17 Jp is clearly visible. 

We compare, in Fig. ^, results of x averaged over dif- 
ferent random realizations, obtained from both the MSW 
theory (14 samples) and QMC simulations (28 samples). 
For a system of 360 spins in 30 FM segments {p — 1/12) 
with AFM couplings (Jaf — O.bJp) between neighbor- 
ing segments, we plot x on both linear scale and log-log 
scale (inset). On the log-log scale, we can easily identify 
the crossover from the FM spin-chain physics to random 
spin-chain physics around = 0.13Ji?- The two meth- 
ods are in good agreement, suggesting that the MSW 
theory can be generalized to a certain class of random 



spin chains. Near the crossover, the MSW method over- 
estimates X, according to the QMC simulations. The 
MSW theory is not exact even in pure FM segments; up 
to 5% error can be observed for a 12-spin FM segment 
near this crossover temperature. In addition, by enforc- 
ing a periodic boundary condition for each FM segment, 
we add an extra FM coupling per segment, which can 
enhance the FM behavior of the system. We require that 
X drop to zero at T 0; the MSW method seems to 
feature a faster drop for such a short chain, compared to 
the QMC simulations. 

The modified spin-wave theory to random spin chains 
has beeru-applied to one-dimensional Heisenberg ferri- 
magnets.im In this case, two constraints corresponding to 
zero-magnetization on each of the two sublattices should 
be applied simultaneously, similar to the multiple zero- 
magnetization constraints on all FM segments in random 
spin chains. The linear combinations of the two chem- 
ical potentials (local fields), in our terminology, have 
the physical interpretations of uniform magnetic field 
and staggered magnetic field, which ensure both mag- 
netization and staggered magnetization be zero. The 
MSW study on ferrimagnetic chaintO revealed that the 
zero-staggered-magnetization constraint keeps number of 
spin waves finite at both high-T and low-T limit, while 
the zero-magnetization constraint does only at the low- 
T limit. We point out that the two chemical potentials 
are equivalent to the two kinds of Lagrangian multipliers 
in the Schwipger-boson mean-field study of the ferrimag- 
netic chains .113 
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